Exact shape of the lowest Landau level in a spin-^ system with uncorrelated disorders 
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Using the path-integral approach developed by Brezin et al. [ Nucl. Phys. B 235, 24 (1984) ], 
we obtain an analytical expression for the density of states of a spin-1 disordered two-dimensional 
electron gas in a strong, perpendicular magnetic field. The density of states of this system illustrates 
the interplay between the Zeeman splitting of Landau levels and the disorder-induced broadening. 
We find that the broadening and the band splitting of the Landau bands are enhanced due to the 
level repulsion from the mixing of two spin orientations caused by the random scatterings. The 
comparison between the spin-| model and the double-layer system is also discussed. 
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I. INTRODUCTION 

There has been remarkable interest in the problem of 
a disordered two-dimensional electron gas in a strong, 
perpendicular magnetic, field ever since the discovery of 
the quantum Hall effcctul in such systems^-Although there 
have been many studies on this systema La, however, only 
a few rigorous results of the density of states (DOS) for 
the disordered electrons are obtained. The difficulty for 
obtaining the rigorous expression of DOS comes from the 
fact that, in the absence of the disorder, the energy spec- 
trum is discrete. Therefore the self-energy of an elec- 
tron is real in any finite order of the perturbation the- 
ory. Hence, to obtain the expected disorder-broadened 
Landau levels ( LL ), one needs to sum over the entire 
diagram expansion. In the strong- field limit, where only 
the lowest LL is considered, the exact DOS of the spin- 
less electrons was first found by Wegnem for the case of a 
Gaussian white noise distribution of the random poten- 
tial. Using a Junctional-integral approach, Brezin, Gross 
and ItzyksonQ were able to re-derive Wegner's result and, 
furthermore, to generalize it to the case in which a non- 
Gaussian, but still zero-range, distribution of impurities 
is assumed. The crucial point for these exact deriva- 
tions is the intrinsic supersymmetry of the problem, by 
which a mapping from the original two-dimensional sys- 
tem onto a zero-dimensional model becomes possible. 
Thus, the summation of the entire diagram expansion 
can be achieved. 

However, in the real world, electrons carry spins. As 
a result, the effects of mixing between LL's with dif- 
ferent spin quantum numbers would become important, 
when the Zeeman splitting is smaller than the width of 
each disorder- broadened LL. Hence, the generalization 
of above results to the spin-i case is of pratical inter- 
est. However, as was pointed out in the last section of 



Ref.EI, the identity ensuring the existence of the super- 
symmetry fails to hold in the two-component case and 
implying that a straightward generalization to the spin-i 
case appears impossible in general. Nevertheless, we find 
that an exact expression of DOS can be reached provided 
that the spin-i electrons are influenced by the random 
potential scattering and the random spin-flip scattering, 
which are both assumed to be the zero-range white-noise 
distribution. We find that the energy difference of the 
splitted bands is larger than the Zeeman splitting en- 
ergy. Moreover, when the Zeeman splitting is turned 
off, the width of the disorder-broadened band is larger 
than, rather than identical to, that of the spinless case. 
Thus an interesting interplay between the Zeeman split- 
ting of Landau levels and the disorder-induced broaden- 
ing can be demonstrated in our spin-i system. In Sec- 
tion II, we introduce the model Hamiltonian and outline 
the functional-integral formulation. In Sec. Ill, the exact 
results of DOS are presented, which show the enhance- 
ment of the splitting and the band width. In Sec. IV, we 
discuss the asymptotic behavior of DOS when Zeeman 
splitting is large. Section V is devoted to the conclusion 
and the comparison between the spin-i model and the 
double-layer system. 

II. THE MODEL AND THE 
FUNCTIONAL-INTEGRAL FORMULATION 

We consider the following Hamiltonian, 

H = H + U 1 (r) + U 2 (r)T 1 , (1) 
Ho = ^(P - eA) 2 - \g m Br z , (2) 

where A(r) = (—By/2,Bx/2) is the vector potential of 
the constant magnetic field, B(r) = Bz, in the symmet- 
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ric gauge, g is the gyromagnetic ratio, /zb is the Bohr 
magneton, t\ and T3 are Pauli matrices. The random 
potential scattering and the random spin-flip scattering 
are denoted by U\(r) and E/2(r)Ti, respectively. The lat- 
ter term can be visualized as a simplified description of 
the scatteijiag by magnetic impurities, or the spin-orbit 
scattering. E£L( A similar model has been considered by 
Wang et alxA to study the effects of Landau level mixing 
)■ 

We first outline briefly the procedure of calculating 
DOS developed by Brezin et alu. Assuming that the 
magnetic field is so strong that the random scatterings 
cannot induce transitions between different Landau lev- 
els, the averaged DOS of the lowest LL can be written 



p(E) = ImTr 

7T 



f 



'E-H + iO' 



(3) 



where Tr denotes the trace operation over the spin indices 
and the bar indicates averaging over all configurations 
of the random potentials. Using the functional-integral 
approach,Q the matrix element of the resolvent after ran- 
dom average becomes 



Tr 



f 



1 E - H + iO 



|r) 



5 -§« 2 M 2 



/ II Du <jDu* a Dv a Dv a 2J ( u l u a + v a v a ) e f 

J ~ U ~ L 



(4) 



a=± <r=± 

with the "action" S given by 



S = 



dzdz*e 2« 2 I Z I 2 ^(e- ag) {u* a u a + v a v a ) 



I 



dzdz* fi 



e -§« 2 N 2 J2 ( 

<7=± 



U„U a 



(5) 



dzdz* f 2 



a=± 



where k 2 = eB/h, e = E — 9 = \9 HbB, and 

z = x + iy denotes the complex coordinates on the plane. 
The holomorphic bosonic ( fermionic ) fields and their ad- 
joints are denoted by u±(z) ( v±(z) ) and u±(z*) ( v±(z*) 
) respectively. The Fourier transform of the distribution 
function Pj(Uj) at a single site, fj(ot), is defined by 



fj (a) = In / dUjPj (U^e-^ , f 3 (0) = 0, 



j = l,2. 
(6) 



In this paper, we restrict ourselves to the case where 
both P\ and P2 are white-noise distributions of the same 
strength, w, i.e., 



Pi{Ui 



exp — 7, — ), 

2w 



P 2 («7 2 ) = exp(-K), (7) 



such that /i(a) = /2(a) = —^a 2 . To examine whether 
the supersymmetry exists in this path-integral represen- 
tation, we introduce the holomorphic supcrficlds and 
their adjoints for each spin orientation, 
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$±(z,6) = u±{z) + -j=k6v±{z), 
$±(z*,6) = u* ± {z*) + -L K v ± (z*) 



(8) 
(9) 



where 9 and 9 are the Grassmann coordinates satisfying 



dOdO 



dOdOO 



dodoe = 0, 



dodeee 



1 

(10) 



Consequently, the " action " S can be expressed in terms of the superfields §±(z, 9) and &±(z* 

d:d : ' ,/YM?e~5« 2 (M 2 +ee) [e($ + $+ + $_$_) - - 
1 ,.2 1 1 2 2tt 



2m 



dzdz* 
dzdz* 



e 2' 



e 2' 



K 
K 2 



d9d9e~^ K ™($+$ + + 



as 



(11) 



Obviously, the first integral in Eq. (|1 1|) is invariant un- 
der superspace rotations, 



5z = lu9, 
59 = ujz, 



Sz* = 6u, 
59 = z*u, 



(12) 
(13) 



which preserves the quadratic term (zz* + 99). However, 
as z and 9 are not treated on the same footing in the sec- 
ond and the third integrals of Eq.(prl), the invariance for 



these integrals is not necessarily promised. For the spin- 
less electrons described by the superfield <I>, the following 
identity 



% I d9d9e~^ 2v& <M> 

K 



n k 



(14) 



ensures the supersymmetry and then the exact DOS of 
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spinless electrons. However, as was pointed out in Ref.B, 
the generalization of Eq. ( |l^) to the two-component case, 
say, $ — » (<&+<&+ + $_$_), is not allowed when n > 1 
as the left-hand side of Eq. (^4|) contains the four fermion 
interactions, v+v+ V—V-, whereas the right-hand side of 
Eq.(|l4|) does not. Hence, the supersymmetry, which en- 
ables one to obtain an exact expression of DOS, does not 
seem to exist in the spin-^ case. However, when the two 
random distributions Pi and P2 are both Gaussian and of 
the same strength, these four fermion interactions cancel 
each other in the sum of the second and third integrals 
in Eq.([ll|). This cancellation can be easily illustrated by 
the following changes of variables for the superfields: 

$1 = -L($ + + 4>_), # 1 = i=($ + + #_), (15) 

$ 2 = -L (*+-$_), # 2 = - $_), (16) 

therefore, a supersymmetric " action " is reached, 



S = — — I dzdz* dOdOe 



l lK , 2 (\z\ 2 +e§) 



x [e($i$i + $2$2) - 3($1$2 + $2$l) 
K 2 J 

X [(#l$l) 2 + (# 2 $ 2 ) 2 ] ■ 



(17) 



In the derivation of Eq.(|l7|), the identity Eq.(|l4"|) is used 
for each $1 and $2 in the intermediate step, in which 
no four fermion term is involved. Because of the dimen- 
sional reduction, d — ► d — 2, due to the existence of the 
supersymmetry in Eq.(^j), the calculations of the aver- 
aged resolvent, Eq.(|), can be reduced to evaluate the 
following zero-dimensional integral,Q 



TW (r 



1 



|r> 



E- H + iO' 
J dtidttdcfadffitylfc + (f>* 2 (f>2)e So 



J d^d^d^d^e^ 
where So is defined by 

2ir 2ir 
So = i-o e(0*0i + $$ 2 ) ~ i—o9(<l>*<h + 4>*2<t>i) 



(18) 



wir 



[(^</> 1 ) 2 + (0^ 2 ) 21 



(19) 



Combining Eqs.(|), © and ©, we have the following 
exact expression for DOS 



p(E) = ^lm^lnZ 



with 



Z 



(20) 



(21) 



It is convenient to rewrite Zq in a different form. We first 
rescale (j>i by <f>\ = \J^:4>i, and subsequently decouple 
the quartic terms in So with the help of gaussian integral 
over a pair of auxiliary variables, Ai and A2. Finally we 
calculate the remaining integral over ^ to obtain 



Zn = 



(*tt) 2 



00 />oo 



7r ( K 2 ) 2 r 2 j v ./ v 

R 1 A ry\ 



d\\d\2- 



Z + iZ 



(e + A 1 )(e + A 2 )-g 2 
(22) 

with the real and the imaginary parts given by 

e~ y2 



7 K 



dy 



x ^ e 
K A 



2 _ 2 / ^ 

ffae* — e + dxe 1 
Jo 



Z = / d ?/ 



-3/ 



2r 2 7 ^y 2 + 25 2 /r 2 



(23) 
(24) 



where the magnetic-field-dependent width T = yfwn 2 /ir : 
v± = v2e/r ± \Jy 2 + 2g 2 /T 2 , and the changes of vari- 
ables, x = (Ai + X 2 )/V2 and j/ = (Ai — A 2 )/\/2, are 
introduced. 



III. RESULTS OF VARIOUS DENSITY OF 
STATES 

The general behavior of DOS is shown in Fig.l. Only 
the part for e > is shown, since p{& is symmetric with 
respect to the band center ( e = ).c3 Figure 1(a) shows, 
as expected, that the splitting of the disordered Lan- 
dau bands increases as g increases. To get more insight 
from our result, we depict in Fig. 1(b) by shifting hori- 
zontally each curve of Fig. 1(a) with the Zeeman energy, 
g = \g hbB. It becomes apparent that, as g increases, 
the shifted curves evolve to a single curve haviipg the simi- 
alar form of DOS in the spinless case, po(E)Ei However, 
two unexpected results emerge: first, the peak values of 
the splitted bands occur at a larger e, rather than exactly 
at e = g as one expects naively; secondly, if there is no 
Zeeman splitting ( i.e., = 0), from Eqs.(22)-(p4[) and 
(pb|) , we obtain 



p{E)\ g=a = V2^ 



2ir 2 T ^ 



1 



(25) 



where To = \J wk 2 /2tt = T/t/2 is the corresponding 
magnetic-field-dependent width in the spinless case. ( 
See Eqs.(46a) and (46b) in Ref.Q. ) Thus the DOS in 
the spin-i system without Zeeman splitting is not sim- 
ply twice as large as that in the spinless case. In our case, 
the peak value at the band center ( e = ) of Eq.([25|) 
is merely y/2 times of that in the spinless case, and the 
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width of p(E) is \/2 times wider. It is because the ran- 
dom scatterings mix the two spin degrees of freedom, as 
indicated by the quartic terms in Eq.(p"l|). Therefore, we 
may expect a level splitting .due to the usual repulsion 
of eigenvalues by these terms,c3 and thus the broadening 
of the Landau bands is enhanced when g = 0, while the 
band splitting becomes larger for g ^ 0. 



- p{E ) ^ -p(E), 



(30) 



and then we have 



FIG. 1(b) 




Fig.l : (a) The averaged DOS for a spin-i system in 
units of the peak value of the spinlcss case, P a(E)\ e= o — 
y/2 K 2 /ir 5 / 2 r, for values of g/T =0.0 ( points ), 0.4 ( dot- 
ted line ), 0.8 ( dashdotted line ), 1.2 ( dashed line ), and 
2.0 ( soild line ). (b) The curves in (a) are shifted hori- 
zontally by the Zeeman splitting energy, g = ^g pbB. 

Moreover, the exact expression of the disorder- 
averaged DOS for each spin component, which is denoted 
by p]{E) ( pi(E) ) for up ( down ) spin, can be obtained 
for the present model. Parallel to the discussion for the 
derivation of the sum of DOS's of two spin components, 
we have 



p(E) = p t (E) + Pl (E). 



(26) 



On the other hand, their difference can be expressed by 
( c.f. EqM) ) 



p(E) = pi(E) - Pl (E) 



K \ 9 1 7 



Consequently, we have 
1 



Pl(E) = ~(p(E)+p(E)), 
Pi {E) = \{p{E)-p{E)). 



(27) 

(28) 
(29) 



By a similar reasoning of the mirror symmetry about the 
band center ( e — ) for p(E) 10 one can easily show that 



p ^E)^ Pi {E), 
Pi(£)™ pi(E). 



(31) 
(32) 



With these relations in mind, we show our results of 
p r (E) and Pi {E) only for e > in Figs.2(a) and (b). We 
find that, as g increases, P[{E) diminishes, while P]{E) 
evolves to a similar form of po(E). The latter can be 
seen more clearly, if we shift horizontally each curve in 
Fig. 2 (a) with g as shown in Fig. 3. The interesting points 
are that : first, the peak value of P](E) ( Pi(E) ) does not 
occur at e = g ( e = ~g ); secondly, P](E) and P[{E) at 
the band center diminish very fast so that a small bump 
of pi(E) ( p-\ (E) ) is developed around e ~ g ( e ~ — g ). 
Again, these results can be explained by the level repul- 
sion from the mixing of two spin degrees of freedom. 




Fig. 2 : The averaged DOS for (a) the spin-up compo- 
nent, Pi(E), and (b) the spin-down component, P[{E), in 
units of the peak value of the spinless case, p {E)\ t=0 = 
V2 K 2 /7r 5 / 2 r, for values of g/T =0.0 ( points ), 0.4 ( dot- 
ted line ), 0.8 ( dashdotted line ), 1.2 ( dashed line ), and 
2.0 ( soild line ). 
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FIG.3 




(e-g)/r 

Fig. 3 : The shifted DOS for the spin-up component, 
Pt(E), in Fig. 2 (a) by the Zeeman splitting energy, g = 
\g p B B. 



IV. ASYMPTOTIC BEHAVIOR OF DENSITY OF 
STATES FOR LARGE ZEEMAN SPLITTING 

In the following, we consider the asymptotic behav- 
ior of DOS in large-g limit ( g/T ^> 1 ) and give the 
analytical descriptions of the features mentioned above. 
Especially, we will focus our attention on the three re- 
gions of energy : (1) near the location of the peak; (2) at 
the band center; (3) toward the band tails. 

Due to the exponential factor, exp(— y 2 ), in the in- 
tegrands of Eqs.(|23]) and (|24|), the characteristic value 
of y, being of order 1, is much smaller than g/T when 
g/T ^> 1. Thus we can neglect y 2 in \Jy 2 + 2g 2 /T, and 
then v± = \plx± , where x± = (e ± g) /T is independent 
of y. In this case, the integrations over y in Eqs.([2^) and 
(E4T) can merely contribute constant factors. Therefore, 
the real part and the imaginary part of Zq become 



2V2Tg 



(/--/+), 



(33) 
(34) 



J ° 2V2Tg 2 
where f± = e~_^f±, and F± — e~ 2x ± Jq* X± dxe x ^ . From 



Eqs.(|(| and (§7|), we have 



p(E) = 
p(E) = 



V2k 2 4y/2f(/ + F--/-F + ) 
7r5/2rl(F_-F + )2 + (/_-/ + )2' 

V2n 2 4^2[f (f+F_ - f_F+) + f(f. 



/+)] 



5/2p 



l(F__ J F + )2 + (/__/ + )s 



(35) 



(36) 



Consequently, the spin-up and spin-down components of 
DOS become 



V2n 2 2^[z_(/ + F_ -I-F+) 



4(f- 



Pi{E) = - 



^ 5/2 r ±(F_-F + y + (f_-f + ) 2 

V2k 2 2V2lx + (f + F_ - f_F + ) + - /+)] 



7r 5 / 2 r 



Hf- 



F+) 2 + (/- 



Uf 



It can be easily shown that the mirror symmetry about 
e = still remains in the large-g limit. 

Let us first consider the behavior of p^(E) near the 
location of its peak value. From Fig. 3, we find that 
the peak of p-\{E) occurs at e = go, which is somewhat 
larger than g, and the deviation of go from g decreases 
as g increases. Hence, p^{E) can be approximated to its 
second-order Taylor polynomial about e = g, and then 
we have 



pm = A{- 



B) 



with 



A = 



B = - 



1 1 



0( 



1 



c 



(S/T) 
V2k 2 



1 



(1 



(3/r) 3 
1 



(9/r) 2 
), 



c 



), 



ir 5 / 2 T y 32(g/T) 2 ' 
where the asymptotic expansion 



dxe x S — (1 
2z v 



(s/r) 4 



2z 2 ' 



(39) 

(40) 
(41) 

(42) 
(43) 



for large z is used for F± in Eq.(|37|). Thus we find 
that the peak of p-\(E) occurs at go/T = g/T + B 



V2n 2 

7r 5 / 2 r 



8(I7r) 



with the peak value p^{E)\ e= g Q = C = 



1 



32(l7rp 



). For <?/r = 2.0, we have g /T = 2.06 



and C = 0.98 , which agrees with the solid line in 
Fig. 2 (a). Moreover, in the limit g — > 00, we have 50 ~~ * <?> 
and /3| (i?) around e — go converges to 



P^E) = 



7r 5 / 2 r n 



d-i)(^) 

7T I 



(44) 



where e' = e — 3. We find that, up to a shift in energy 
by g, our result is identical to po{E) near its band center 
(c.f. Eq.(48) in Ref.LJ). Therefore, in the large-g limit, 
the DOS of each spin component in our spin-^ system 
seems to behave as the shifted po(E). However, we will 
show below that this is not the case. 

Now let us turn to the region at the band center. For 
e = 0, i.e., x± = ±g/T, we have /+ = /_ = er 2 ^/ T )\ 

•v%/r 



and F- 



+ 

follows that 



-F- 



Pl(E)U=o = 



r 

Jo 



T/_ 9 
' 1 



dxe x . By Eq.(43), it 



7r 5 / 2 r 2 v r 



,2 e -(g/r ) 



(45) 



which is only one-half of the -value of po(E)\e=±g f° r 
g/T Q > 1 (c.f. Eq.(47) in Ref.fl). It indicates that the 
behavior of pi (E) ( and pi (E) , due to the mirror symme- 
try ) far from the peak is not identical to that of Pq(E) 
in the band tail. 

We can show the difference in the asymptotic be- 
havior between po(E) and Pi(E) ( pi(E) ) by consid- 
ering the regions around the band tails directly. For 
(e — q)/r > q/r 3> 1, namely, x + 3> X— 1, from 
Eqs.(p7|) and (|43|), we obtain 



Pl(E) £* 



7r 5 /2r 2 v r 



(^) 2 ( : 



2 ^--(e7r ) 2 



(46) 



Due to the mirror symmetry about e = 0, the behavior 
of P]{E) in the limit e' — > — oo ( or equivalently, e — > — oo 
) can be inferred by that of p\ (E) in the limit e' — > oo. 
Therefore, by Eqs.(38) and ((43|), we have 



K 7T.I 



_^5)2 e -(e7r ) 2 



(47) 



which diminishes faster than that in the opposite side of 
the band tail. Hence, the asymptotic behavior in both 
limits, d — ► ±oo, does not reduce to that of the shifted 
po(E) in the band tails ( c.f. Eq.(47) of RefE ), as in- 
dicated in the discussion for e = 0. Moreover, we find 
that the asymptotic behavior of pi(E) ( pi(E) ) is not 



symmetric with respect to the position of its peak. Thus 
Pi (E) and pi (E) do not evolve to the identical form of 
the shifted p n (E) in all regions of energy, although they 
do so near the location of their peaks. 



V. CONCLUSION AND DISCUSSION 



In summary, we extend the exact solution obtained by 
Brezin et al. for the DOS of spinless electrons at the 
lowest Landau level with a short-range disorder to the 
cases of a spin-i system. In this generalization, we ob- 
tain an analytical expression for DOS of spin-i system, 
illustrating the interplay between the Zeeman splitting of 
Landau levels and the disorder-induced broadening. We 
also find that the broadening and the band splitting of 
the Landau bands are enhanced due to the level repul- 
sion from the mixing of two spin degrees of freedom by 
the random scatterings. 

We notice that a similar expression of DOS is recently 
reached by Shahbazyan and RaikhEHl in the context of the 
double-layer system with a tunneling coupling constant 
t. Alternatively, we can compare these two systems ex- 
plicitly by calculating their DOS using the path-integral 
approach, rather than summing up the entire diagrams 
in Rcf.lH3. After introducing the holomorphic superficlds 
and their adjoints ( i = 1, 2 ) for electrons in each 
layer, the " action " of the double-layer system, Syy, be- 
comes ( c.f. Eq.(O) ) 



St 



2ni 



J dzdz*d6d6e~? K2(lzl2+0S} [e(*i*i + #2*2) - t(tfi#2 + #2*1) 



dzdz*gi 
dzdz* gi 



K 2 J 



(48) 



e-W\A^JL I d9d ee 



where g\{a) and ^2(0) are the corresponding Fourier 
transforms of the distribution functions for the random 
potentials in each layer. For the white-noise distributions 
of the same strength, by Eq. (|4|) , So reduces to the same 
form as our supersymmetric " action " in Eq.(p"^), and 
then leads to a similar expression of DOS. Thus, by the 
change of variables in Eqs.(|l5|) and (|l6|), our model can 
be related to the double-layer one. However, the equiva- 
lence between these two models is not held in general. If 
the random distributions are not white-noise type of the 
same strength, the cancellation of the four fermion inter- 
actions will fail, and then there is no supersymmetry for 
the spin-i case. However, for the double-layer system, a 
supersymmetric " action " for Su can still be obtained ( 
because there is no four fermion interaction in Su ) , and 
an exact expression of DOS can be formulated for any 
zero-range random distributions by a-straightward gen- 
eralization of the work of Brezin et aln That is, once the 



disorder is included, these two systems are not exactly 
the same and, moreover, for the spin-i model, it is quite 
limited to have an exact expression of DOS. 
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If we set e 



-e, then v± 



we have -> , Z\ 



and. 



by Eqs.(||) 



and 



^(Zl/Zf). Thus, after rewriting Eq.@ as p(E) = 
£°-(Zl/Z£)/[l + (Zl/Z$f], the symmetry with re- 
spect to e — > — e becomes obvious. 
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